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Discussions  of  energy  dissipation  during  friction  processes  have  captured  the  attention  of  engineers 
and  scientists  for  over  300  years.  Why  then  do  we  know  so  little  about  either  dissipation  or  friction 
processes?  A  simple  answer  is  that  we  cannot  see  what  is  taking  place  at  the  interface  during  sliding. 
Recently,  however,  devices  such  as  the  atomic  force  microscope  have  been  used  to  perform  friction 
measurements,  characterize  contact  conditions,  and  even  describe  the  “worn  surface.”  Following 
these  and  other  experimental  developments,  friction  modeling  at  the  atomic  level — particularly 
molecular  dynamics  (MD)  simulations — has  brought  scientists  a  step  closer  to  “seeing”  what  takes 
place  during  sliding  contact.  With  these  investigations  have  come  some  answers  and  new  questions 
about  the  modes  and  mechanisms  of  energy  dissipation  at  the  sliding  interface.  This  article  will 
review  recent  theoretical  and  experimental  studies  of  friction  processes  at  the  atomic  scale. 
Theoretical  treatments  range  from  simple,  analytical  models  of  two-dimensional,  coupled 
ball-spring  systems  at  0  K,  to  more  complex  MD  simulations  of  three-dimensional  arrays  of 
hydrogen-  and  hydrocarbon-terminated  surfaces  at  finite  temperatures.  Results  are  presented  for  the 
simplest  yet  most  practical  cases  of  sliding  contact:  sliding  without  wear.  Sliding  without  friction  is 
seen  in  weakly  interacting  systems.  Simple  models  can  easily  explain  the  energetics  of  such  friction 
processes,  but  MD  studies  are  needed  to  explore  the  dynamics  (excitation  modes,  energy 
pathways,...)  of  thermally  excited  atoms  interacting  in  three-dimensional  fields.  These  studies 
provide  the  first  atomic-scale  models  for  anisotropic  friction  and  boundary  lubrication.  Friction 
forces  at  atomic  interfaces  must  ultimately  be  measured  at  the  macroscopic  level;  these 
measurements,  which  depend  on  the  mechanical  properties  of  the  measuring  system,  are  discussed. 
Two  rather  unique  experimental  studies  of  friction  are  also  reviewed.  The  first  employs  a  “surface 
force  apparatus”  to  measure  adhesion  and  friction  between  surfactant  monolayers.  The  correlation 
of  adhesion  hysteresis  and  friction  provides  a  new  mechanism  of  friction;  moreover,  the 
interpretation  for  the  effect — hysteresis  from  entanglement  of  the  molecular  chains  during  a  phase 
transformation — implies  that  the  dynamics  are  taking  place  at  an  accessible  time  scale  (seconds  to 
minutes).  The  second  study  extends  the  time  domain  at  which  friction  can  be  measured  to  the 
nanosecond  scale.  A  quartz-crystal  oscillator  is  used  to  monitor  the  viscosity  of  monolayer  liquids 
and  solids  against  solid  surfaces.  Interfaces  slip  angstroms  in  nanoseconds.  Modelers  have 
suggested  a  variety  of  mechanisms  for  this  atomic-scale  friction  process,  from  defect-mediated 
sliding  to  electron  drag  effects.  The  article  ends  by  identifying  the  vast,  barely  charted  time-space 
domain  (micro-to-pico  time  and  length  scales)  in  which  experiments  are  needed  to  further 
understand  the  dynamic  aspects  of  friction  processes. 


I.  BACKGROUND 

Friction  can  now  be  studied  at  the  atomic  scale,  thanks  to 
developments  in  the  past  decade  of  a  variety  of  experimental 
techniques.1  The  most  well  known,  often  referred  to  as  proxi¬ 
mal  probes,  have  evolved  from  scanning  tunneling  micros¬ 
copy  (STM);2  they  include  atomic  force  microscopy  (AFM)3 
and  its  sliding  companion  friction  force  microscopy, 
(FFM).4  6  These  probes  allow  friction  to  be  studied  with 
atomic  resolution  in  all  three  dimensions.  Another  proximal 
probe,  generically  known  as  a  surface  force  apparatus  (SFA), 
affords  atomic  resolution  only  in  the  vertical  direction,  but 
allows  direct  measurement  and/or  control  of  micrometer¬ 
sized  areas  of  contact  in  the  lateral  direction. 7-14  A  very  re- 


a*E-mail:  singer@chem.nrl.navy.mil 


cent  technique,  based  on  the  quartz  crystal  microbalance 
(QCM),  permits  sliding  friction  processes  to  be  studied  at  the 
angstrom  level  and  at  time  scales  in  the  nanosecond 

15-17 

range. 

Although  friction  processes  may  originate  at  the  sliding 
interface,  the  measurement  of  friction  is  usually  performed 
by  macroscopic  devices — springs,  levers,  dashpots,  etc. — 
often  located  far  from  the  interface.  In  order  to  link  measured 
frictional  forces  with  theory,  it  has  been  necessary  to  exam¬ 
ine  the  influence  of  mechanical  parameters,  such  as  stiffness, 
on  friction  measurements,  not  an  unfamiliar  problem  to 
tribologists.18'19  Proximal  probes,  however,  are  sensitive  to 
mechanical  properties  of  the  probe  as  close  as  the  first  atomic 
layer  of  the  tip  and  as  far  away  as  the  compliant  lever 
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The  opening  of  experimental  studies  of  friction  at  the  na¬ 
nometer  and  nanosecond  scale  has  attracted  theorists 
equipped  to  model  physical  and  chemical  processes  at  these 
scales.  Surface  scientists  are  now  using  sophisticated  solid- 
state  potentials  to  calculate  mechanical  interactions  between 
surfaces.22  Friction  force  calculations  are  performed  either 
analytically  or  by  molecular  dynamics  (MD)  simulations.23'24 
MD  simulation  affords  the  added  opportunity  of  using  video 
animation  to  study  friction  processes.  For  the  first  time,  sci¬ 
entists  can  “see  what  is  happening”  at  the  otherwise  buried 
sliding  interface.  As  you  will  read  shortly,  they  see  atomic 
and  molecular  excitation  modes  (vibrations,  bending,  and  ro¬ 
tations),  electron-hole  excitations,  density  waves,  and  mo¬ 
lecular  interdigitation,  to  name  a  few.  Once  the  modes  are 
identified,  physicists  and  chemists  can  address  perhaps  the 
most  fundamental  but  least  understood  aspect  of  friction:  en¬ 
ergy  dissipation  processes.  While  it  has  been  known  for 
centuries25'26  that  most  frictional  energy  dissipates  as  heat, 
neither  the  macroscopic  nor  microscopic  mechanisms  of  en¬ 
ergy  dissipation  have  been  fully  explained.  Adhesion  and 
deformation  contributions  to  energy  dissipation  in  friction 
have  recently  been  discussed  by  Johnson.27  Flere  I  review 
some  recent  theoretical  and  experimental  studies  of  atomic- 
scale  friction  behavior  and  modes  of  energy  dissipation. 

Before  launching  into  the  studies,  it  is  useful  to  review  the 
thermodynamic  criteria  used  to  study  energy  dissipation. 
Clearly,  if  friction  processes  generate  heat,  they  are  irrevers¬ 
ible  and  cannot  be  treated  by  classical  thermodynamics.  If, 
however,  each  step  in  a  sliding  interaction  is  executed  with 
infinite  slowness  and  with  the  two  couples  always  in  equilib¬ 
rium  (never  an  unbalanced  force),  then  the  process  may  be 
considered  reversible.  In  such  a  quasistatic  process,  sliding 
can  be  achieved  with  zero  friction.  Real  systems  can  ap¬ 
proximate  such  reversible,  adiabatic  processes  so  long  as  the 
rate  at  which  each  step  is  taken  is  much  slower  than  the 
relaxation  time  of  the  system.28  However,  any  instability  in 
the  mechanical  system  that  leaves  unbalanced  forces  will  re¬ 
sult  in  an  irreversible  process  in  which  energy  is  lost  (or, 
more  precisely,  unrecoverable).  In  a  mechanical  system, 
where  force  F  acts  over  a  distance  r,  the  energy  lost  over  a 
cycle  is  given  by 

A U=  <j>  F  dr.  (1) 

In  an  atomic  sliding  calculation,  a  cycle  can  be  a  translation 
across  some  periodic  distance  of  the  lattice,  e.g.,  one  atomic 
spacing.  In  an  experiment,  one  cycle  can  be  a  single  pass 
over  a  surface,  including  the  making  and  breaking  of  the 
contact. 

This  article  is  presented  in  the  following  sequence.  The 
second  section  deals  with  theoretical  approaches  used  to  ex¬ 
amine  friction  processes  and  friction  measurements  at  the 
atomic  scale.  Conditions  that  can  give  zero  friction  and  finite 
friction  in  simple  systems  are  presented.  Two  molecular  dy¬ 
namics  simulations  of  more  realistic  yet  “wearless”  friction 
studies  are  described  and  modes  of  energy  dissipation  are 
identified.  The  third  section  presents  two  experimental  ap¬ 
proaches  that  have  succeeded  in  identifying  microscopic  fric¬ 
tion  process;  in  both  cases,  energy  losses  are  identified  with 


hysteresis  losses  in  the  system.  The  fourth  section  summa¬ 
rizes  our  understanding  of  interfacial  friction  processes  and 
energy  dissipation  mechanisms,  and  the  fifth  section  consid¬ 
ers  ways  that  atomistic  approaches  can  be  used  to  solve  prac¬ 
tical  problems  in  tribology. 

II.  THEORETICAL  APPROACHES 
A.  Frictionless  sliding 

Two  analytical  studies  that  identify  conditions  for  zero 
friction  and  the  transition  to  finite  friction  are  reviewed.  The 
first  investigates  the  potential  seen  by  a  moving  atom,  while 
the  second  examines  the  force  experienced  by  the  tip  of  a 
FFM.  Although  both  present  surface  interactions  at  0  K,  they 
establish  base  line  criteria  for  zero  friction  systems  and  mea¬ 
surements  of  zero  friction  by  which  more  general  calcula¬ 
tions  can  be  evaluated. 

McClelland29  30  describes  the  sliding  friction  behavior  of 
a  simple  two-dimensional  couple  consisting  of  two  sub¬ 
strates:  the  stationary  upper  substrate  has  an  atom  BQ  at¬ 
tached  to  a  spring  while  the  moving  lower  substrate  consists 
of  equally  spaced,  rigid  atoms  (see  Fig.  1).  Because  of  the 
periodicity  of  the  lower  layer,  the  entire  sliding  behavior  can 
be  analyzed  by  following  the  motion  of  atom  B0  across  one 
atomic  spacing.  According  to  this  independent  oscillator  (IO) 
model,  atom  B0  experiences  forces  exerted  by  a  spring  above 
and  the  atoms  below;  the  forces  are  derived  from  potentials 
VBB  and  V/AB,  respectively.  As  the  lower  substrate  moves,  the 
combined  interaction  potential  Vs  changes.  The  changing  po¬ 
tential  and  the  atomic  trajectory  of  atom  B{)  are  depicted  in 
the  five  sketches  running  left  to  right  across  the  figure.  The 
upper  sketches  represent  a  “strong”  interaction  and  the 
lower,  a  “weak”  interaction.  The  atom’s  position  on  a  spring 
is  shown  as  an  open  circle  and  its  position  in  the  potential  V s 
by  a  solid  circle.  In  a  “strong”  potential,  the  atom  is  initially 
repelled  to  the  right;  then  beyond  half  an  atomic  spacing,  it 
“snaps”  back  to  the  left,  the  atomic  equivalent  of  “stick- 
slip.”  The  snap  back,  or  “plucking”  motion,  can  be  under¬ 
stood  in  terms  of  the  evolving  shape  of  potential  Vs.  In  an 
adiabatic  process,  an  atom  must  always  sit  at  a  minimum  in 
the  potential  well.  As  can  be  seen,  however,  the  position  of 
the  atom  at  (d)  is  only  a  local  minimum,  which  disappears  by 
position  (e).  At  the  moment  the  local  metastable  minimum 
“flattens  out,”  there  are  unbalanced  forces  on  atom  B0;  to 
regain  an  equilibrium  state,  the  atom  falls  to  the  position  of 
the  stable  minimum.  Since  this  transition  does  not  occur  un¬ 
der  equilibrium  conditions,  the  process  is  irreversible  and  the 
energy  lost  in  the  fall  cannot  be  reused  to  assist  the  sliding 
process.  Stated  in  more  physical  terms,  the  strain  energy  put 
into  stretching  the  spring  at  B0  is  not  recovered  locally;  in¬ 
stead,  it  is  converted  into  vibrational  motion  which  dissipates 
into  the  substrate  (as  heat). 

This  instability  can  be  avoided  by  using  a  “weaker”  in¬ 
teraction  potential.  As  shown  in  the  lower  sequence  of 
sketches  in  Fig.  1,  the  weaker  potential  does  not  develop  a 
local  minimum.  The  atom  moves  smoothly  through  a  repul¬ 
sive  then  an  attractive  force  field,  first  being  repelled  by,  then 
pushing,  the  lower  substrate.  Since  the  system  remains  in 
equilibrium  throughout  the  cycle,  no  energy  is  dissipated  and 
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Fig.  1 .  Two  representations  of  the  motion  of  atom  B0  in  the  independent  oscillator  model:  atom-on-spring  and  atom-in-potential.  Top  and  bottom  rows  depict 
strong  and  weak  interfacial  interactions,  respectively.  The  left-most  diagrams  display  the  relevant  potentials  (see  text);  subsequent  panels  illustrate  the  response 
of  B0  to  progressive  sliding  of  the  lower  layer  of  atoms.  The  location  of  atom  B0  is  represented  by  a  black  dot  in  the  combined  potential  Vs  plotted  below  each 
atom-on-spring  diagram  (Refs.  29  and  30). 


the  friction  force  is  zero.  McClelland  then  gives  more  precise 
criteria  for  stability  and  shows  that  qualitatively  similar  fric¬ 
tion  behavior  occurs  with  other  more  complex  sliding 
couples. 

Tomanek  efa/.31-33  also  describe  conditions  associated 
with  frictionless  sliding,  emphasizing  the  mechanical  proper 
ties  of  the  apparatus  as  well  as  the  strength  and  shape  of  the 
interatomic  potentials.  They  present  two  idealized  models  for 
a  FFM  tip  interacting  with  an  atomic  surface,  only  one  of 


which  is  presented  here.  Called  a  “realistic-friction  micro¬ 
scope”  by  the  authors,  the  model  accounts  for  the  elasticity 
of  the  FFM  as  well  as  the  external  load  Fext  and  the  surface 
interaction  potential. 

Figure  2  shows  a  FFM,  with  a  horizontal  spring  that  pulls 
the  spring-tip  assembly  along  the  interface  potential  of  the 
substrate.  The  tip  experiences  the  combined  force  of  the  in¬ 
terface  potential  and  Fcn.  As  the  tip  moves  across  the  sur- 


Fig.  2.  Model  of  a  FFM.  External  suspension  M  is  guided  along  a  horizontal 
surface  in  the  x  direction.  The  load  Fext  is  kept  constant  along  the  trajectory 
shown  by  arrows.  The  spring-tip  assembly  is  elastically  coupled  to  the  sus¬ 
pension  M  in  the  horizontal  direction  by  a  spring  of  constant  c.  The  friction 
force  Ft  is  related  to  the  elongation  xt—xM  of  the  horizontal  spring  from  its 
equilibrium  value  (Refs.  31  and  32). 


Fig.  3.  The  friction  force  Ft  as  a  function  of  the  FFM  position  xM  for  a  hard 
and  soft  spring  and  the  average  friction  force  (Ff)  for  a  soft  spring  (Refs.  31 
and  32). 
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Fig.  5.  Friction  vs  temperature  data  in  normalized  units  for  the  case  of 
strong  f 6 1 / 6(j  =  1 .0)  and  weak  l6|/e0=0. 1 1  interfacial  interactions  at  a  sliding 
speed,  V=20  m/s.  The  solid  line  represents  a  fit  to  the  data  using  a  thermal 
activation  model  (Ref.  35). 


Fig.  4.  Contour  plot  of  the  average  friction  force  ( F t- )  per  atom  as  a  function 
of  the  load  Fext  and  the  horizontal  spring  constant  c  (for  monoatomic  Pd  tip 
on  graphite)  (Refs.  31  and  32). 


face,  the  horizontal  spring  elongates  by  an  amount,  xM—xt, 
which  depends  on  the  stiffness  c  of  the  horizontal  spring.  For 
a  “hard"  spring  (c>ccrit),  both  the  tip  trajectory  and  the  F f 
curve  are  single-valued  functions  of  xM\  an  example  of  the 
latter  curve  is  labeled  “hard”  in  Fig.  3.  Since  the  positive 
and  negative  excursions  of  F f  are  the  same  over  a  cycle,  the 
average  friction  force  (F  f)  is  zero.  When  c  falls  below  ccrit, 
both  the  tip  trajectory  and  the  F y  curve  are  triple-valued 
functions  of  xM;  an  example  of  the  latter  is  labeled  “soft”  in 
Fig.  3.  In  this  “soft”  spring  case,  the  tip  snaps  forward  at  the 
point  of  “instability,”  as  in  the  “strong”  spring  case  in  the 
IO  model  above,  giving  the  asymmetric  F f  versus  xM  curve 
shown  by  the  curve  with  arrows.  Flere,  the  average  friction 
force,  given  by  the  dashed  line,  is  nonzero.  The  energy  dis¬ 
sipated  by  the  collapse  of  the  elongated  spring,  calculated 
according  to  Eq.  (1),  is  represented  by  the  shaded  area  under 
a  portion  of  the  F f  versus  xM  curve. 

Figure  4  gives  the  average  friction  for  a  range  of  values  of 
horizontal  spring  constant  c  and  external  force  Fext.  Two 
results  are  apparent.  First,  the  ability  to  measure  a  zero  fric¬ 
tion  process  depends  on  the  value  of  c.  Second,  the  friction 
coefficient  /jL=(Fj)/Fext  for  atomic  contacts  is  not  indepen¬ 
dent  of  load.  Further  discussion  of  friction  versus  load  be¬ 
havior  for  atomic  sliding  is  given  elsewhere.31,32 

These  simplified  models  of  friction  between  atomic 
couples  provide  analytical  criteria  for  transitions  to  zero  fric¬ 
tion,  thus  zero  energy-dissipation  conditions  at  T =0  K.  Zero 
friction  requires  “weak”  interaction  forces  (low  atomic- 
scaled  corrugations)  to  achieve  adiabatic  motion  and  “hard” 
horizontal  springs  to  measure  the  effect.  An  atomic  “stick- 
slip”  phenomena  is  predicted  when  these  criteria  are  not  met. 
The  models  are,  in  fact,  consistent  with  atomic  scale  FFM 
measurements.4 

B.  MD  simulations  of  monolayer  films 

MD  simulations  of  friction  behavior  go  beyond  the  simple 
analytical  models.  They  use  more  realistic  interaction  poten¬ 


tials  that  exhibit  anisotropy  in  two  or  three  dimensions  and 
they  account  for  effects  of  temperature  and  sliding  velocity. 
Moreover,  video  animations  of  the  simulations  allow  us  to 
visualize  the  trajectories  of  the  atoms  at  and  near  the  inter¬ 
face.  MD  thereby  gives  us  our  “first  look”  at  what  happens 
at  the  buried  interface  during  sliding  contacts.  Early  ex¬ 
amples  of  MD  studies  of  frictional  contact  between  solid 
surfaces  depicted  wear  and  transfer  of  material  in  the  sliding 
contact.24,34  More  recently,  two  studies  of  “wearless”  fric¬ 
tion  have  been  reported;  “wearless”  friction  means  that  no 
atoms  are  lost  from  or  transferred  to  either  member  of  the 
couple.  In  both  studies,  the  solid  surfaces  were  terminated 
with  a  monolayer  of  a  simple  hydrocarbon  or  hydrogen.  This 
section  presents  these  studies,  which  begin  to  address  the 
role  of  surface  films  in  friction  processes. 

1.  Monolayers  of  alkane  chains 

McClelland  and  Glosli30,35  have  performed  MD  simula¬ 
tions  of  friction  between  two  monolayers  of  alkane  chains. 
The  chains,  six  carbon  atoms  long,  are  initially  ordered  in  a 
herringbone  pattern.  Chain  bonds  are  allowed  to  bend  and 
twist,  but  not  stretch.  Carbon  and  hydrogen  atoms  on  each 
chain  interact  with  atoms  of  other  chains  by  a  Lennard-Jones 
potential  of  interaction  strength  e0.  The  interaction  strength 
at  the  interface  el  is  adjustable  in  order  to  study  its  influence 
on  friction  behavior.  The  temperature  of  the  layers  is  held 
constant  by  means  of  a  heat  bath;  energy  losses  in  the  chains 
can  thus  be  followed  as  heat  losses  in  the  layers.  MD  calcu¬ 
lations  are  performed  over  a  temperature  range  of  0 <T <300 
K  and  sliding  velocities  up  to  t/(l=204  m/s.  The  friction  force 
is  calculated  as  the  shear  stress  t  averaged  over  several 
cycles  of  sliding. 

Figure  5  shows  two  sets  of  friction  versus  temperature 
data  for  the  case  of  strong  (el/eQ=  1.0)  and  weak  (e1/e0=0.1) 
interactions  at  a  sliding  speed  of  V=20  m/s.  For  both  cases, 
the  friction  force  exhibits  different  behaviors  in  the  three 
(low,  medium,  and  high)  temperature  regimes;  in  addition, 
for  the  case  of  weak  interactions,  the  friction  vanishes  at  low 
and  high  temperatures.  In  order  to  interpret  the  friction  be¬ 
haviors  and  energy  dissipation  processes,  McClelland  and 
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Fig.  6.  Shear  stress  (left)  and  heat  flow  (right)  vs  sliding  displacement  for 
strong  (elleQ=  1.0)  and  weak  (  £|/e0=0. 1 )  interfacial  interactions  at  7’ =20  K. 
The  “zero”  heat  flow  line,  in  lower  right  box,  was  added  by  present  author 
(Ref.  35). 


Glosli  relied  on  video  animations  of  the  molecular  dynamics 
and  calculations  of  both  shear  stress  and  heat  flow  versus 
displacement  curves. 

At  low  temperatures,  the  trajectories  of  the  alkane  chains 
are  very  much  like  that  of  the  atom-on-a-spring  in  the  IO 
model  discussed  earlier.  For  strong  interactions,  the  ends  of 
the  chains  are  strained  until  the  local  minimum  in  the  poten¬ 
tial  disappears  and  the  ends  are  released  abruptly.  The  chains 
oscillate  with  a  pivoted-hinge  motion,  retaining  their  herring¬ 
bone  pattern.  At  this  temperature,  chain  oscillations  cause  the 
CHX  groups  on  the  backbone  to  vibrate  but  do  not  induce 
torsional  (twisting)  motion  of  the  chains  themselves,  i.e.,  vi¬ 
bration  and  twisting  modes  are  decoupled  at  low  tempera¬ 
tures. 

Friction  force  versus  distance  curves  exhibit  atomic  stick- 
slip,  as  seen  in  the  upper  left  of  Fig.  6.  At  slip,  strain  energy 
is  released  abruptly  (25  ps),  which  caused  the  heat  flow 
shown  in  the  upper  right  of  Fig.  6.  Because  the  energy  dis¬ 
sipation  is  associated  with  a  plucking  instability,  the  friction 
is  expected  to  be  independent  of  velocity;  this  was  confirmed 
by  MD  calculations  at  several  velocities  for  T=20  K.  Pluck¬ 
ing  motion  and  energy  dissipation  occur  until  the  interfacial 
interaction  ratio  falls  below  €1/e0=0.4.  Below  this  value,  the 
friction  force  versus  distance  curves  vary  smoothly  and  sym¬ 
metrically  about  F/= 0  and  there  is  no  measurable  heat  flow 
(see  lower  left  and  right  parts  of  Figs.  6,  respectively).  Like 
the  weakly  interacting  atom  in  the  IO  model,  the  weakly 
coupled  alkane  chains  exhibit  harmonic  motion  as  the  two 
substrates  slide  past  each  other. 

As  the  temperature  increases,  the  friction  in  both  interac¬ 
tion  ranges  increases,  reaching  maximum  values  around 
T=80  K.  At  this  temperature,  video  animations  depict  very 
complex  dissipation  modes.  The  strain  energy  released  by  the 
bent  chain  now  couples  into  three  different  modes:  chain 
oscillations,  CH3  group  vibrations  and,  for  the  first  time,  the 
torsional  (twisting)  modes.  Torsional  modes,  which  are  “fro¬ 
zen”  at  lower  temperatures,  become  excited  at  intermediate 
temperatures;  this  process  is  sometimes  referred  to  as  “rota¬ 
tional  melting.”  With  all  three  modes  now  anharmonically 


coupled,  excitations  damp  out  quickly  as  the  energy  is  trans¬ 
ferred  to  lattice  vibrations  in  the  substrate.  Hence,  the  in¬ 
crease  in  friction  at  intermediate  temperatures  can  be  attrib¬ 
uted  to  the  excess  energy  associated  with  the  unfrozen 
torsional  modes. 

As  the  temperature  increases  above  the  rotational  melting 
temperature,  a  third  mechanism  comes  into  play:  molecules 
vibrate  so  actively  in  all  directions  that  an  increasing  per¬ 
centage  of  the  chains  can  hop  and  slide  over  the  opposite 
surface  without  introducing  strain.  This  reduces  the  net  fric¬ 
tional  force  in  both  weak  and  strong  interaction  cases;  in 
fact,  in  the  weak  interaction  system,  the  friction  force  goes  to 
zero  at  highest  temperatures. 

Glosli  and  McClelland  then  demonstrate  that  the  friction 
behavior  at  high  temperatures  is  consistent  with  classical  be¬ 
havior  of  polymeric  films.  MD  calculations  of  strongly  inter¬ 
acting  couples  showed  increasing  friction  (sublinear)  with 
sliding  velocity  (not  shown  here)  and  decreasing  friction 
with  temperature  (Fig.  5).  The  latter  curve  was  fitted  to  the 
Eyring  thermal  activation  model  that  Briscoe  and  Evans36 
used  to  describe  the  shear  behavior  of  thin  polymeric  films. 
The  straight-line  fit,  drawn  in  Fig.  5,  required  only  one  ad¬ 
justable  parameter,  the  activation  energy  Q\  it  was  deter¬ 
mined  to  be  Q  =70  K.  The  shear  strength  (at  T= 0  K)  derived 
from  Q  was  r=32  MPa,  about  the  same  value  obtained  ex¬ 
trapolating  to  T= 0  K  Briscoe  and  Evans’  experimental  data 
for  stearic  acid  (C18)  and  behenic  acid  (C22). 

In  summary,  MD  simulations  of  friction  behavior  between 
monolayers  of  alkane  chains  show  that  the  simple  harmonic 
versus  plucking  friction  behavior  applies  to  more  compli¬ 
cated  systems;  in  addition,  the  simulations  enumerate  the 
modes  of  dissipation  that  arise  with  increasing  molecular 
complexity  and  increasing  thermal  activation. 

2.  H-  and  (H+CxHy)-terminated  diamond  (111)  surfaces 

Harrison,  White,  Colton,  and  Brenner37  39  have  simulated 
the  friction  behavior  of  H-  and  (H+CtHv)-terminated  dia¬ 
mond  (111)  surfaces  placed  in  sliding  contact.  The  forces  are 
derived  from  an  empirical  hydrocarbon  potential  capable  of 
modeling  chemical  reactions  in  diamond  and  graphite  lattices 
as  well  as  small  hydrocarbon  molecules.40  Two  diamond 
(111)  surfaces,  each  terminated  in  a  (1 X 1)  pattern,  are  placed 
in  twin  (mirror-image)  contact;  the  distance  of  separation 
depends  on  the  repulsive  interaction  potential  and  the  load. 
Sliding  is  performed  in  two  directions:  the  [112]  direction 
and  the  [110]  direction.  In  the  [112]  direction,  opposing  hy¬ 
drogen  atoms  can  make  “head-on”  contact  at  very  high  loads 
because  their  velocity  vectors  lie  in  a  common  plane  perpen¬ 
dicular  to  the  surface;  this  will  be  called  the  “aligned"  di¬ 
rection.  In  the  [110]  direction,  the  hydrogen  atoms  can  only 
pass  adjacent  to  each  other  (across  each  others’  diagonals) 
because  their  velocity  vectors  do  not  lie  in  a  common  per¬ 
pendicular  plane.  The  lattice  temperature  is  set  at  300  K, 
unless  otherwise  stated,  and  sliding  velocities  are  50  or  100 
m/s.  Normal  loads  are  varied  up  to  0.8  nN/atom,  correspond¬ 
ing  to  mean  pressures  up  to  20  GPa,  and  friction  coefficients 
are  averaged  over  a  unit  cell. 

Their  first  study  examines  two  H-terminated  diamond 
(111)  surfaces.37  Along  the  [112]  direction,  the  friction  coef- 
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Fig.  7.  Average  friction  coefficient  as  a  function  of  normal  load  for  (H 
+2A)-terminated  surfaces  sliding  in  the  [112]  direction  at  v=l  A/ps  and  at 
T= 300  K.  Curves  are  for  X=hydrogen,  methyl,  ethyl,  and  n-propyl- 
terminations,  respectively.  From  Ref.  39. 

ficient  begins  near  zero  for  lowest  loads,  increases  nearly 
linearly  with  load  up  to  0.6  nN/atom,  then  levels  out  at 
fi=0.4  (see  Fig.  7). 

Video  animation  sequences  of  sliding  along  the  “aligned” 
[112]  direction  show  different  interaction  mechanisms  at  low 
and  high  loads.  At  lowest  loads,  opposing  hydrogen  atoms 
first  repel  each  other  backwards.  Then  as  strain  develops, 
they  pivot  sideways  and  revolve  past  each  other  at  closest 
approach,  finally  pushing  each  other  forward.  With  these  tra¬ 
jectories,  the  net  frictional  force,  averaged  over  a  unit  cell, 
nearly  cancels  and  almost  no  energy  is  dissipated.  At  higher 
loads,  atomic-level  stick-slip  occurs:  instead  of  gently  re¬ 
volving  by  each  other,  opposing  hydrogen  atoms  momen¬ 
tarily  become  “stuck”  in  a  (mutually)  repulsive  potential 
well,  then  suddenly  “slip”  and  revolve  around  one  another. 
The  pivoting-then-revolving  motion,  assisted  by  thermal  mo¬ 
tion,  excites  both  vibrational  and  bending  modes  in  the  C-H 
bonds.  These  excitations  are  passed  on  to  the  lattice  as  vibra¬ 
tions  (phonons)  then  heat.  In  this  way,  potential  (strain)  en¬ 
ergy  developed  at  high  loads  is  transformed  into  kinetic  en¬ 
ergy  (mechanical  excitations),  leading  to  nonzero  values  of 
the  friction  coefficient. 

Friction  coefficient  versus  load  data  for  sliding  along  the 
[110]  direction  (not  shown  here)  are  about  an  order-of- 
magnitude  lower  than  along  the  [112]  direction.  Animations 
show  that  hydrogen  atoms,  instead  of  meeting  along  aligned 
directions,  “zigzag”  through  channels  of  potential  minima 
that  run  between  adjacent  rows  of  hydrogen  atoms  on  the 
opposing  surfaces.  Since  the  opposing  hydrogen  atoms  do 
not  encounter  each  other  directly,  strain  levels  are  lower  and 
thus  the  frictional  forces  are  lower.  Hence,  different  sliding 
directions  on  identical  surfaces  lead  to  anisotropy  in  both 
atomic  trajectories  and  friction  coefficients.  Perhaps  these 
trajectories  can  account  for  the  well-know  frictional  anisot¬ 
ropy  of  single  crystals.41,42 


The  friction  coefficient  of  the  H-terminated  surface  shows 
a  temperature  dependence  that  is  qualitatively  similar  to  that 
found  with  alkane  chains.  At  a  fixed  pressure  of  3  GPa,  the 
friction  coefficient  at  T= 0  K  is  high  (/x=0.4),  but  drops  as 
the  temperature  rises:  to  /x=0.25  at  70  K,  then  to  /x=0. 1 5  at 
300  K.  The  friction  coefficient  is  larger  at  low  temperatures 
because  opposing  hydrogen  atoms  cannot  rotate  out  of  the 
way  of  aligned-trajectory  repulsive  wells  without  the  help  of 
thermal  motion. 

Harrison  et  air*  have  also  examined  the  friction  behavior 
of  the  same  sliding  configuration  of  H-terminated  diamond 
but  with  two  methyl  groups  substituted  for  two  hydrogen 
atoms  on  one  surface.  Along  the  [112]  direction,  the  friction 
coefficient  versus  load  data  rise  quicker  than  on  the  fully 
H-terminated  surface,  but  reach  lower  steady-state  values  (fi 
=0.35  versus  /r=0.4)  (see  Fig.  7).  These  differences  can  be 
attributed  to  the  large  volume  occupied  by  a  methyl  group. 
Instead  of  easily  “revolving”  around  the  hydrogen  atoms 
during  low-load,  aligned-trajectory  collisions,  the  methyl 
group  gets  stuck,  then  “slips”  with  a  ratcheting  motion.  The 
ratcheting  motion,  analogous  to  the  motion  of  a  “turnstile,” 
rotates  the  methyl  group  alternately  clockwise  then  counter¬ 
clockwise  120°  around  the  C-C  bond.  Turnstile  motion,  ac¬ 
companied  by  C-H  bond  excitations,  is  responsible  for 
higher  friction  at  low  loads.  At  high  loads,  the  size  of  the 
methyl  groups  keeps  the  two  surfaces  further  apart  than  com¬ 
parably  loaded  H-terminated  pairs.  Harrison43  speculates  that 
the  “flattening”  out  of  the  friction  coefficient  versus  load 
data  may  be  due  either  to  screening  of  the  interaction  poten¬ 
tial  or  to  constraining  the  excitation  modes. 

Sliding  the  methyl-substituted  surface  along  the  [110]  di¬ 
rection  (not  shown)  produces  the  same,  strong  load- 
dependent  friction  coefficients  found  in  the  [112]  direction 
(see  Fig.  7).  Remarkably,  the  substitution  of  two  CH3  mol¬ 
ecules  for  two  hydrogen  atoms  produces  an  order  of  magni¬ 
tude  greater  friction  coefficient.  Why?  Unlike  the  terminal 
hydrogen  atoms,  which  can  “zigzag”  freely  through  the  ad¬ 
jacent  hydrogen  atom  channels,  the  larger  methyl  groups  ex¬ 
hibit  “turnstile”  rotations  like  those  found  in  [112]  sliding. 
The  rotations  are  accompanied  by  “zigzag”  motion,  which 
becomes  more  pronounced  as  the  load  increases.  The  in¬ 
creased  energy  expended  by  the  larger  molecules  in  this  turn¬ 
stile  and  zigzag  trajectory  is  responsible  for  the  increase  in 
friction  coefficient  with  load. 

In  their  latest  study,  Harrison  et  al.i9  have  substituted  two 
ethyl  and  two  « -propyl  groups  for  two  hydrogen  atoms;  fric¬ 
tion  coefficient  versus  load  data  are  seen  in  the  two  right 
panels  of  Fig.  7.  The  larger,  more  flexible  hydrocarbon 
groups  reduce  friction  at  high  loads  by  a  factor  of  1.5  to  2 
compared  with  the  fully  H-terminated  diamond  slider. 
Center-of-mass  trajectories  of  the  CH3  portion  of  the  ethyl 
groups,  plotted  on  potential  energy  contour  maps  of  a 
H-terminated  diamond  (111)  surface,  give  insight  into  how 
the  motion  of  an  ethyl  molecule  affects  the  friction  coeffi¬ 
cient.  At  low  loads  (not  shown  here),  the  ethyl  molecule 
bends  over,  lies  down,  and  is  dragged  almost  straight  across 
the  repulsive  potentials,  like  the  trajectory  a  chain  would 
have  if  one  end  were  tied  to  the  upper  surface.  At  high  loads, 
however,  the  ethyl  molecule  uses  its  flexibility  and  length  to 
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Fig.  8.  Center-of-mass  trajectories  of  the  CH3  tail  of  an  ethyl  group  (dark 
solid  line)  plotted  on  a  potential  energy  contour  map  of  a  H-terminated, 
diamond  (111)  surface.  The  two  solid  triangles  indicate  the  starting  points 
for  two  CH3  tails.  Sliding  is  along  the  Y  direction.  The  two  dashed  lines 
represent  the  trajectories  of  the  attachment  point  out  the  ethyl  groups  during 
sliding.  Average  normal  force  is  0.42  nN/atom.  See  Ref.  39  for  details. 

“snake”  (detour)  around  high  potential  energy  barriers.  This 
trajectory,  depicted  by  solid  lines  in  Fig.  8,  expends  less 
energy  and  produces  a  lower  friction  coefficient  than  the 
pivot-rotation  or  turnstile  modes  found  with  smaller  mol¬ 
ecules. 

In  summary,  Harrison  el  al.  have  identified  several 
mechanisms  which  may  account  specifically  for  the  friction 
behavior  and  energy  dissipation  of  hydrogen-  and 
hydrocarbon- terminated  diamond  surfaces,  and,  more  gener¬ 


ally,  for  boundary  film  lubrication.  They  have  shown  that,  at 
300  K,  hydrogen-  and  hydrocarbon-terminations  follow  dif¬ 
ferent  trajectories  when  sliding  along  “hard”  (aligned- 
trajectory)  and  “soft”  (channels)  directions  of  H-terminated 
diamond,  and  thereby  explain  the  strong  frictional  anisotropy 
of  H-terminated  diamond  pairs  along  selected  low-friction 
channels.  They  have  cataloged  numerous  excitation  modes 
(rotations,  turnstiles,  ...)  by  which  frictional  energy  is  dissi¬ 
pated.  In  addition,  they  have  shown  that,  at  high  contact 
stresses,  larger  hydrocarbon  groups  reduce  friction  even  fur¬ 
ther  because  of  size  and  steric  effects. 

III.  EXPERIMENTAL  APPROACHES 
A.  Friction  and  adhesion  hysteresis 

Israelachvili  and  co-workers44  have  recently  discovered  a 
new  relationship  between  adhesion  and  friction,  based  on 
experimental  studies  of  surfactant  monolayers.  Experiments 
are  performed  with  a  SFA,  in  which  both  adhesion  and  fric¬ 
tion  are  measured.  Adhesion  behavior  is  examined  in  contact 
radius  versus  load  curves  during  loading -unloading  cycles 
and  in  pull-off  force  measurements;  friction  behavior,  in  uni¬ 
directional  and  reciprocating  sliding.  The  surfactant  mono- 
layers  studied  exhibit  one  of  three  phases:  solidlike,  amor¬ 
phous,  or  liquidlike.  The  amorphous  state  is  a  phase  in 
between  the  solidlike  and  liquidlike  state.  Moreover,  the 
phases  of  each  of  the  layers  could  be  changed  by  varying  the 
atmosphere,  temperature,  velocity  or  related  parameters. 

The  main  conclusion  of  the  study  is  that  the  friction  force 
does  not  correlate  with  the  adhesion  force  (or  adhesion  en- 


Inert  Air  Decane  Vapor 


Fig.  9.  (A)  and  (B)  Friction  traces  of  two  CaABS  monolayers  at  25  °C  exposed  to  inert  air  and  to  air  saturated  with  decane  vapor.  (C)  and  (D)  Adhesion 
energies  on  loading,  unloading  and  pull-off  measured  under  the  same  conditions  as  the  upper  friction  traces  (Ref.  44). 
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Fig.  10.  A  schematic  “friction  phase  diagram”  representing  the  trends  ob¬ 
served  in  the  friction  forces  of  five  different  surfactant  monolayer  types 
studied.  The  curve  also  correlates  with  adhesion  hysteresis  of  the  monolay¬ 
ers  but  not  with  the  adhesion  per  se  (Ref.  44). 


ci'gv  y),  but  rather  with  hysteresis  in  the  adhesion  force.  As 
an  example.  Fig.  9  shows  friction  and  adhesion  measure¬ 
ments  for  loading  and  unloading  two  calcium  alkylbenzene- 
sulfonate  (CaABS)  monolayers  at  25  °C.  The  curves  to  the 
left  (A  and  C)  represent  behavior  of  a  liquidlike  monolayer: 
quite  a  low  friction  coefficient  and  some  hysteresis  in  the 
adhesion  during  the  loading -unloading  cycle.  The  curves  to 
the  right  (B  and  D)  show  that  after  exposing  the  surfaces  to 
decane  vapor,  the  already  low  friction  coefficient  decreases 
even  more  and  the  adhesion  hysteresis  disappears.  Previous 
studies45'46  had  shown  that  when  hydrocarbon  vapors  con¬ 
densed  onto  CaABS,  the  molecules  penetrated  the  outer 
chain  regions  and  fluidized  the  surface.  Thus,  it  was  hypoth¬ 
esized  that  if  the  monolayer  were  made  more  liquidlike,  the 
friction  and  adhesion  hysteresis  would  be  reduced. 

Many  other  correlations  between  friction  and  adhesion 
hysteresis  and  the  phase  of  surfactant  monolayers  have  been 
observed.  Both  friction  and  adhesion  hysteresis  increase 
when  solidlike  monolayers  or  liquidlike  monolayers  are 
made  amorphous.  Conversely,  when  amorphous  monolayers 
are  made  more  fluidlike  or  solidlike,  both  the  friction  and 
adhesion  hysteresis  decrease. 

Observed  trends  in  friction  and  adhesion  hysteresis  be¬ 
havior  are  summarized  in  the  schematic  “friction  phase  dia¬ 
gram  curve”  shown  in  Fig.  10.  Maximum  values  of  friction 
and  adhesion  hysteresis  ( but  not  adhesion  values)  are  found 
around  a  “chain-melting”  temperature  Tm;  this  is  the  tem¬ 
perature  at  which  the  monolayer  is  inbetween  the  solidlike 
and  liquidlike  state,  i.e.,  the  amorphous  state.  Lower  values 
of  friction  and  adhesion  hysteresis  are  found  at  temperatures 
above  or  below  Tm.  It  is  seen  that  the  amorphous  state  rep¬ 
resents  the  highest  friction  and  highest  adhesion  hysteresis. 
Factors  that  can  change  the  phase  state  of  monolayers,  such 
as  vapors,  speed,  etc.,  can  effectively  shift  the  curve  in  di¬ 
rections  indicated  by  the  arrows  in  Fig.  10. 

Israelachvili  et  a/.44  give  a  physical  basis  for  this  behav¬ 
ior.  Adhesion  hysteresis  is  the  “irreversible”  part  of  the  ad¬ 
hesion  energy,  and  is  related  to  the  energy  dissipated  during 


the  loading  and  unloading  process.  A  likely  molecular  origin 
of  adhesion  hysteresis  is  the  extent  of  interpenetration  and 
subsequent  ease  of  disentanglement  of  the  molecules  across 
an  interface.  If  there  is  little  interpenetration,  as  with  solid¬ 
like  layers,  the  friction  is  smooth  and  no  additional  energy  is 
expended  separating  surfaces.  If  there  is  significant  interpen¬ 
etration,  as  with  liquidlike  layers,  but  also  ease  of  disen¬ 
tanglement  on  separation,  the  friction  is  again  low  and  little 
extra  energy  is  expended  separating  surfaces.  A  thermody¬ 
namic  description  of  the  liquidlike  case  would  conclude  that 
the  time  to  separate  the  chains  is  slower  than  the  relaxation 
time  of  the  molecules  and,  therefore,  that  separation  approxi¬ 
mates  a  reversible  process.  In  both  the  above  cases,  the  sys¬ 
tems  are  physically  in  similar  states  going  into  and  coming 
out  of  contact.  By  contrast,  with  amorphous  layers,  there  is 
significant  chain  interpenetration  but  separation  occurs  faster 
than  the  molecular  relaxation  time.  Thus,  amorphous  layers 
will  be  in  different  states  going  into  and  coming  out  of  con¬ 
tact,  and  consequently,  more  energy  will  be  expended  to 
separate  them  than  would  be  needed  for  the  two  other 
phases. 

This  new  mechanism  of  friction  in  boundary  lubrication 
can  operate  over  a  wide  time  scale,  the  time  it  takes  a  par¬ 
ticular  molecule  to  adapt  its  trajectory  to  the  lowest  possible 
interaction  potential,  relative  to  the  time-dependence  of  the 
potential.  In  the  case  of  the  MD  studies  of  n  -propyl- 
terminated  diamond,39  the  low  friction  was  achieved  at  high 
loads  because  the  molecules  could  “follow”  the  lowest  en¬ 
ergy  contours  in  picoseconds.  The  same  was  true  for  the 
liquidlike  monolayers,  whose  chains  could  follow  the  mini¬ 
mum  force  trajectories  in  milliseconds.44  However,  on  the 
same  time  scale,  amorphous  chains  could  not  follow  mini¬ 
mum  force  trajectories,  and  energy  differences  were  seen  in 
both  adhesion  hysteresis  and  friction. 

B.  Friction  and  slip  of  monolayer  films 

Krim  et  al. 15-17  are  studying  the  frictional  forces  for  sol¬ 
idlike  and  liquidlike  films  adsorbed  on  conducting  (metal) 
and  insulating  (oxidized  metal)  substrates.  Their  approach  is 
quite  novel.  Films  of  gases  such,  as  Kr  and  Ar  or  C2H4  and 
C2H6,  are  condensed  onto  surfaces  to  thicknesses  up  to  sev¬ 
eral  monolayers.  The  thickness  of  the  film  is  determined  in  a 
straight-forward  manner  with  a  QCM.  In  addition,  the  QCM 
monitors  subangstrom  shifts  in  the  vibrational  amplitude 
caused  by  gas  adsorption.  These  shifts  are  due  to  frictional 
shear  forces  between  the  condensing  film  and  the  oscillating 
surface.  Krim  et  al.15  have  shown  that  the  “slip  time”  rof  a 
monolayer  film  can  be  determined  with  subnanosecond  ac¬ 
curacy  from  these  shifts.  Note  that  the  time  and  length 
scales,  nanoseconds  and  angstroms,  makes  these  experiments 
unique  in  the  field  of  tribology.  Since  slip  is  fundamentally 
an  energy-dissipative  process,  the  technique  allows  energy 
dissipation  to  be  measured  and  the  mechanisms  of  energy 
dissipation  to  be  studied. 

Experiments  with  rare  gas  atoms  have  shown  that  (1)  the 
slip  times  for  monolayers  physisorbed  on  smooth  gold  sur¬ 
faces  are  on  the  order  of  nanoseconds  and  (2)  solidlike  films 
exhibit  longer  slip  times  than  liquidlike  films.  These  results 
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are  consistent  with  a  frictional  force  proportional  to  the  slid¬ 
ing  velocity,  indicating  a  viscous  friction  mechanism. 

Three  models  have  been  proposed  to  account  for  energy 
dissipation  of  the  sliding  monolayers.  The  first  two  postulate 
that  phonons  carry  away  the  energy.  Sokoloff,  using  an  ana¬ 
lytical  model  of  friction,47,48  suggests  that  defects  between 
the  incommensurate  monolayer- solid  interface  can  account 
for  the  slip  times.  Robbins  et  al.4<>  have  used  molecular  dy¬ 
namic  simulations  to  determine  the  viscous  coupling  be¬ 
tween  a  driven  substrate  and  an  adsorbed  monolayer  film. 
Requiring  no  arbitrary  parameters,  the  model  gives  excellent 
agreement  with  many  of  the  experimental  observations:  it 
gives  the  correct  magnitude  of  slip  time  r,  a  friction  force 
proportional  to  r  for  physisorbed  films;  and  less  slip  in  liquid 
films  than  in  solid  films.  They  have  also  shown  that  their 
results  agree  with  a  simple  analytic  model,  closely  related  to 
that  of  Sokoloff.  The  slip  time  is  directly  related  to  equilib¬ 
rium  properties  of  the  film:  it  is  proportional  to  the  lifetime 
of  longitudinal  phonons  and  inversely  proportional  to  the 
square  of  the  density  oscillations  induced  by  the  substrate. 

A  third  model,  presented  by  Persson  et  al.,50  postulates 
energy  dissipation  by  electron-hole  scattering.  It  assumes 
that  electrons  in  the  metal  substrate  experience  a  drag  force 
equal  in  magnitude  to  the  force  required  to  slide  the  adsorbed 
film.  This  force  is  estimated  from  measurements  of  the 
change  in  resistivity  of  metal  films  as  a  function  of  gas  cov¬ 
erage.  Calculated  slip  times  are  in  good  agreement  with  ex¬ 
perimental  values  for  adsorbed  rare  gases  and  hydrocarbon 
molecules.  The  model  also  predicts  different  slip  times  for 
C2H4  and  C2H6  adsorbate  films,  but  only  if  electronic  contri¬ 
butions  are  present,  e.g.,  with  metals  but  not  insulators. 
Krim’s  group51  has  recently  tested  this  prediction  by  measur¬ 
ing  slip  times  for  C2H4  and  C2H6  on  silver  and  on  oxygen- 
coated  silver.  They  find  different  slip  times  on  Ag,  but  the 
same  slip  time  on  oxygen-coated  silver,  consistent  with  the 
predictions.  Thus,  based  on  these  rather  unique  studies  of 
friction,  it  appears  that  both  phonon  and  electron  mecha¬ 
nisms  contribute  to  energy  dissipation. 

IV.  SUMMARY  AND  DISCUSSION 

From  studies  just  reviewed  and  others  in  the  literature,  our 
understanding  of  interfacial  friction  processes  and  energy 
dissipation  mechanisms  can  be  summarized  as  follows.  Low 
friction,  including  zero  friction,  can  be  achieved  at  low 
loads,  with  weak  interface  interactions  and  with  “small”  at¬ 
oms  at  the  interface.  The  mechanical  principle  that  explains 
this  behavior  follows  from  the  simple,  one-dimensional,  IO 
model:  the  strain  energy  transmitted  by  interfacial  atoms  dur¬ 
ing  the  first  half  of  the  cycle  is  returned  to  them  during  the 
second  half-cycle.  This  behavior  is  also  observed  in  the  more 
realistic,  three-dimensional  MD  simulations.  The  third  di¬ 
mension  itself  contributes  an  additional  friction  reduction 
channel;  it  provides  the  interfacial  atoms  an  extra  degree  of 
freedom — to  move  out  of  their  common  plane — to  escape 
“stick”  events  along  aligned  directions.  For  example, 
H-terminated  atoms  can  rotate  around  each  other  or  “zigzag” 
along  potential  minima  channels;  these  trajectories  are  not 
available  to  atoms  described  in  two-dimensional  models.  In 
principle,  many  “low  friction  trajectories”  can  be  found 


along  selected  directions  in  real  crystals  having  anisotropic 
interaction  potentials  (corrugations);  these  possibilities  are 
treated  more  quantitatively  by  Hirano  and  Shinjo.52  In  prac¬ 
tice,  however,  too  soft  a  spring  in  the  measuring  device  can 
lead  to  friction  force  instabilities,  resulting  in  measured  fric¬ 
tion  forces  that  are  higher  than  expected. 

Friction  is  increased  by  many  factors.  Strong  interfacial 
interactions  (corrugations),  according  to  the  simple  IO 
model,  give  a  finite  static  friction  force,  then  stick-slip  mo¬ 
tion  between  atoms.  Three-dimensional  potentials  also  show 
atomic-scale  stick-slip  processes,  but  produce  modes,  e.g., 
turnstile  motion,  that  are  more  complex  than  the  equivalent 
two-dimensional  plucking  mode.  Atomic-scale  stick-slip  pro¬ 
cesses  have  been  seen  in  FFM  measurements,  but  at  present 
the  modes  responsible  have  not  been  identified  because  of 
the  relatively  slow  response  time  of  friction  devices.53  An- 
harmonic  coupling  of  excited  modes  establishes  multiple 
pathways  for  energy  dissipation,  thereby  increasing  friction 
coefficients.  An  example  is  the  MD  simulation  of  alkanes  at 
intermediate  temperatures,  where  torsional  modes  become 
allowed,  providing  a  new  pathway  for  energy  dissipation. 
Other  excitation  modes  that  enhance  friction  and  dissipate 
energy  are  density  oscillations,  defects  at  interfaces  and 
electron-hole  coupling  in  metals.  Commensurate  lattices 
have  been  shown  to  increase  friction  forces  by  many  orders 
of  magnitude.47  Friction  force  is  expected  to  increase  with 
increasing  external  force  (load).  However,  as  Zhong  and 
Tomanek33  have  shown,  surface  interactions  can  be  per¬ 
turbed  over  a  selected  load  range,  thereby  lowering  the  fric¬ 
tion  coefficient  as  the  load  increases. 

Temperature  can  influence  friction  behavior  in  several 
ways.  Thermal  activation  of  an  energy-dissipating  mode,  like 
rotational  melting,  increases  friction.  In  contrast,  thermal  ac¬ 
tivation  can  lower  potential  barriers  and  increase  tunneling, 
thereby  reducing  friction.  At  high  temperatures,  thermal  ef¬ 
fects  can  dominate  friction  processes,  giving  liquidlike  (vis¬ 
cous)  friction  instead  of  solid  friction  (finite  static  friction  at 
all  velocities).  The  MD  simulation  of  alkane  chain  friction 
showed  this  transition  from  solid  to  viscous  behavior  with 
increasing  temperature. 

The  size  and  shape  of  molecules  can  also  influence  fric¬ 
tion  behavior.  Small  atoms  or  molecules  may  follow  low- 
friction  trajectories  whereas  larger  atoms  or  molecules  may 
not  “fit”  into  the  same  channels,  resulting  in  higher  friction. 
The  H-  versus  CH3-terminated  diamond  along  the  [110]  di¬ 
rection  is  such  an  example.  By  contrast,  larger  molecules  can 
reduce  friction  more  effectively  than  smaller  molecules  at 
higher  loads  if  they  have  sufficient  flexibility  to  spread 
across  the  surface;  an  example  is  the  high  load  behavior  of 
the  hydrocarbon-terminated  diamond  surfaces.  This  steric  ac¬ 
commodation,  however,  can  increase  friction  if  the  molecule 
cannot  follow  the  minimum-energy  trajectory  as  fast  as  the 
surfaces  move  apart.  Chain  entanglement  between  amor¬ 
phous  films  observed  in  SFA  studies  is  an  example  of  steric 
effects  causing  increased  energy  dissipation. 

A  new  concept  of  friction  behavior  has  been  demonstrated 
by  Israelachvili  et  al.  that  energy  dissipation  is  maximum 
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when  the  time  (and  length)  scales  of  contact  (externally  con¬ 
trolled)  match  the  intrinsic  time  and  length  scales  of  molecu¬ 
lar  interactions.  This  concept  is  consistent  with  thermody¬ 
namic  considerations  of  two  bodies  coming  into  contact.  As 
mentioned  in  the  introduction,  the  degree  to  which  the  con¬ 
tact  process  approximates  a  reversible,  quasistatic  process 
depends  on  the  rate  at  which  each  step  is  taken  compared  to 
the  relaxation  time  of  the  system.  Put  in  terms  of  the  driven 
oscillator  analog,  deviations  from  equilibrium  and  energy 
dissipation  are  maximum  when  time  and  length  scales  of  the 
system  and  driver  are  matched.  However,  classical  thermo¬ 
dynamics  is  not  really  suitable  for  the  treatment  of  contacting 
surfaces  let  alone  sliding  surfaces.  Even  in  the  mildest  con¬ 
tact  circumstances,  in  which  the  two  bodies  retain  their  iden¬ 
tity  after  separating,  equilibrium  was  never  achieved;  at  best, 
the  two  bodies  reached  metastable  equilibrium.  A  more  pre¬ 
cise  description  of  the  thermodynamics  of  contacting  sur¬ 
faces  is  needed. 

Finally,  these  studies  can  give  us  some  new  insights  into 
the  role  of  surface  films  in  friction  processes.  Generalizing 
the  studies  of  Harrison  et  al.,  we  see  that  films  in  which 
“small”  atoms  chemisorb  one-to-one  with  the  substrate  lat¬ 
tice  might  provide  the  screening  needed  to  prevent  interface 
“welding”  and  give  low-friction  trajectories  along  weak  cor¬ 
rugation  channels.  Films  made  of  larger  molecules,  with 
lower  compressibility,  might  reduce  friction  at  high  loads  by 
providing  atomic  screening  as  well  as  steric  accommodation. 


V.  RECOMMENDATIONS 

This  review  was  meant  to  introduce  tribologists  to  some 
of  the  more  recent  investigations  of  energy  dissipation  pro¬ 
cesses  in  interfacial  friction.  Many  recommendations  for  fu¬ 
ture  research  in  atomic-scale  tribology  can  be  made  based  on 
these  preliminary  investigations.  For  example,  the  remarks  in 
the  previous  paragraph  suggest  two  approaches  for  modeling 
boundary  lubricant  films:  (1)  gas  or  solute/additive  interac¬ 
tions  with  surfaces  can  be  modeled  with  small,  single-atom 
terminations  and  (2)  “run-in”  boundary  films  can  be  mod¬ 
eled  by  more  complex  molecular  attachments. 

Many  of  the  ideas  discussed  in  this  paper  were  presented 
at  a  two  week  long  NATO  ASI  meeting  held  in  Braunlage, 
Germany  in  August  1991.  Considerable  time  was  spent  dis¬ 
cussing  “future  issues”  and  suggested  approaches  for  re¬ 
search  in  this  field;  these  have  been  published  in  the  Epi¬ 
logue  to  the  conference  proceedings.54  Here  I  summarize 
only  a  few  of  these  items. 

(1)  How  can  atomistic  modeling  continue  to  make  an  impact 
on  understanding  friction?  on  understanding  lubrication? 

(2)  Can  algorithms  (e.g.,  hybrid  methods)  be  developed  to 
simulate  friction  processes  at  time  and  length  scales 
longer  than  can  be  treated  in  molecular  dynamics  calcu¬ 
lations  alone,  e.g.,  that  extend  computational  simulations 
from  the  nm/fs  scale  to  the  pm!  ps  scale. 

(3)  Can  lubricants  be  tailored  to  take  advantage  of  the  dy¬ 
namic  properties  of  certain  fluids,  e.g.,  the  “chemical 
hysteresis”  of  monolayer  films  discussed  by  Israelachvili 
et  al. 


(4)  In  practical  machines,  sliding  is  sustained  on  surface 
films — whether  organic  lubricants,  oxides,  or  other  solid 
films.  Can  molecular  simulations  help  us  to  understand 
the  chemistry  of  film  formation,  the  mechanical  proper¬ 
ties  of  these  films  and  how  the  films  break  down? 

One  of  the  newest  issues  that  tribology  must  deal  with  is 
the  concept  of  matching  time  and  lengths  scales  in  friction 
studies.  As  we  saw  above,  energy  dissipation  hence  friction 
is  intimately  linked  to  time  and  length  scales.  Moreover,  the 
atomic/molecular  modes  of  interfacial  interactions  operate  at 
time  and  length  scales  far  shorter  than  traditional  tribology 
measurements.  In  the  section  of  the  Epilogue54  entitled  New 
ways  of  probing  friction  processes,  we  asked  “How  can  we 
use  the  power  of  microscopic  modeling  to  gain  new  insights 
into  macroscopic  friction  processes  and,  ultimately,  to  solve 
technological  problems?”  Goddard55  suggests  that  this  can 
be  done  by  progressing  along  the  “chain-linked”  ladder,  il¬ 
lustrated  in  Fig.  11,  from  quantum-level  studies  to  engineer¬ 
ing  design.  His  “hierarchy  of  modeling  tribological  behav¬ 
ior”  unites  atomistic  models,  which  operate  in  very  short 
length-time  scales,  with  engineering  models,  which  describe 
tribological  behavior  in  length-time  scales  observable  by 
more  traditional  measuring  equipment.  This  approach  “...al¬ 
lows  consideration  of  larger  systems  with  longer  time  scales, 
albeit  with  a  loss  of  detailed  atomic-level  information.  At 
each  level,  the  precise  parameters  (including  chemistry  and 
thermochemistry)  of  the  deeper  level  get  lumped  into  those 
of  the  next.  The  overlap  between  each  level  is  used  to  estab¬ 
lish  these  connections.  This  hierarchy  allows  motion  up  and 
down  as  new  experiments  and  theory  lead  to  new  under¬ 
standing  of  the  higher  levels,  and  new  problems  demand  new 
information  from  the  lower  levels.” 

But  where  are  the  experimental  approaches  for  investigat¬ 
ing  the  “lower  (short  scale)  levels?”  As  illustrated  in  Fig.  11, 
most  “friction  machines,”  including  the  proximal  probe  de¬ 
vices,  are  operated  at  long  time  scales.  An  abbreviated  search 
of  recent  literature  produced  only  three  “tribology”  tests  and 
a  fourth  proximal  probe  method  that  come  close  to  investi¬ 
gating  friction  behavior  at  short  time  and  length  scales.  La¬ 
beled  (l)-(4)  in  the  Fig.  11,  they  are  described  here  briefly. 

(1)  Bair  et  al.56  have  used  fast  IR  detectors  to  measure  flash 
temperatures  during  high  speed  frictional  contacts  of  as¬ 
perities  of  length  10  pm  and  greater,  with  time  resolu¬ 
tion  of  about  20  ps. 

(2)  Spikes  et  al.  have  developed  real  time  optical  techniques 
for  investigating  the  physical  behavior  of  EHL  films 
down  to  5  nm  thick57  and  chemical  processes  occurring 
in  contacts  10  pm  wide  by  80  nm  thick.'5 

(3)  Krim  et  al.15'16  have  used  the  quartz  crystal  microbal¬ 
ance  experiments  (described  earlier)  for  probing  atomic 
vibrations  amplitudes  between  0.1  and  10  nm  and  time 
scales  from  10~12  to  10~8  s. 

(4)  Hamers  and  Markert59  have  shown  that  STM  images  are 
sensitive  to  the  recombination  of  photoexcited  carriers 
whose  lifetimes  are  in  the  picosecond  range. 

Clearly,  innovative  experimental  approaches  for  measur¬ 
ing  friction  processes  at  short  and  intermediate  time-length 
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Fig.  1 1 .  Time  and  length  scales  of  present-day  models  and  experiments  in  Tribology.  From  Ref.  54. 


scales  are  needed  to  assist  the  modelers  who  are  already 
there.  Tribologists  should  seek  out  physicists  and  chemists 
working  in  these  time-space  domains  and  form  collabora¬ 
tions  to  carry  out  tribology-oriented  experiments.  As  1  have 
tried  to  emphasize  in  this  article,  much  of  the  progress  in  the 
field  comes  when  experiments  and  theories  can  overlap  on 
the  same  time-length  scales. 

Finally,  there  is  another  time  scale  to  consider,  and  that  is 
“the  question  of  time  to  translating  this  (fundamental) 
knowledge  into  engineering  practice.”  Professor  Dowson 
(Leeds  University,  UK)  told  us  how  this  can  be  accomplished 
in  the  final  part  of  the  Epilogue.60 

“We  have  heard  quite  a  lot  about  the  subject  of  friction 
from  the  atomic  scale  up  to  the  macroscopic  scale.... 

There  are  other  aspects  of  scale  I  think  we  should  reflect 
on.  One  is  the  question  of  time  in  terms  of  translating  this 
knowledge  into  engineering  practice...  The  time  scales  are 
generally  enormous.  ...I  think  you  should  take  note  that  it  is 
going  to  be  about  A  DECADE  IF  NOT  A  GENERATION 
before  that  impact  will  be  seen  fully  in  engineering.  This 
makes  it  all  the  more  important  for  ...  physicists,  chemists 
and  engineers  (to  meet),  so  that  the  engineers  can  absorb  by 
osmosis  the  concepts  that  you  are  revealing....  It  is  important 
that  these  ideas  feed  into  our  consciousness  so  that  we  apply 
them  sensibly  in  future  developments. 

...I  hope  equally  that  scientists  will  not  be  too  frustrated 
by  the  fact  we  take  10  to  20  years  to  incorporate  their  bright 
new  ideas  in  designing  better  skis  or  whatever  it  might  be. 
Let  us  perpetuate  this  interaction  between  groups  of  people 
who  all  have  one  objective  and  that  is  to  understand  the  laws 
of  physics  in  order  to  apply  them  as  effectively  as  we  possi¬ 


bly  can  for  the  good  of  society  through  the  manufacture  of 
reliable,  efficient  and  sensible  products.” 
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